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Abstract. We have used UV (FOCA) and FIR (IRAS) images of six nearby late type galaxies to study the radial 
variation of the UV extinction (deduced from the FIR/UV ratio). We compare the UV extinction gradient with 
other extinction indicators (Balmer decrement) and search for a relation between the extinction, the metallicity 
and the gas surface density among our galaxies. We detect in our small sample a clear relation between extinction 
and metallicity. These observed relations are used to calibrate an empirical recipe useful for extinction correction 
in the UV, visible and near-infrared images of late type galaxies. 

Key words. Galaxies: spiral - Ultraviolet: galaxies - ISM:dust, extinction 



1. Introduction 

The physical properties of late type galaxies, and in partic- 
ular the ones that are determined from the stellar emission 
can be studied only after an accurate dust extinction cor- 
rection. It is thus crucial to have a sound understanding 
of the dust extinction in late type galaxies before studying 
nearby as well as high redshift galaxies (especially since 
the latter are observed in the optical in UV-rest frame, 
where the extinction is particularly severe). 

The definition of an accurate model useful for cor- 
recting the UV, optical and near infra-red data for the 
dust extinction in disc galaxies is particularly difficult 
because of several uncertainties: i) little is known about 
the relative geometrical distribution of the dust and stars 
of various spectral types and ages at small and large 
scales. On small scales, young stars spend a finite time 
within the star f orming regions before migrating out of 
them (e.g. [chariot & Fall 2000; Panuzzo ct al. 2003). On 
large scales, the properties and distributions of HII re- 
gions in spiral discs vary with the morphological type 
(Kcnnicutt et al. 19^). ii) the extinction law (variation 
of the extinction with wavelength) presents large varia- 
tions for different lines of sight in the Galaxy, and different 
shapes in the Galaxy and in the Magellanic Clouds, in p ar- 
ticular in the UV domain (e.g. lSavaee fc MathislllO?^ . It 
is thus unclear which extinction law should be applied, iii) 
Despite recent significant improvement in c haracterising 
the extinction curve (e.g. iDesert et aLlll990|) . the nature 



of the dust (composition, size distribution) is still poorly 
known. 

There are however several observational ways to di- 
rectly estimate the dust extinction in well defined spec- 
tral ranges: using the Balmer lines by comparing the ob- 
served to the expected decrement ijLeaueux et al.l Il98ll: 
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The Balmer decrement gives a good measure of the dif- 
ferential extinction between the lines outside of the ionized 
region. Its use however requires an accurate determination 
of the underlying Balmer absorption. Ideally, one needs 
high quality spectra and good resolution for this purpose. 
In resolved galaxies, this method can only be applied to 
HII regions. Given their peculiar na ture, the effects of age 
and geometrical distribution (e.g. ICharlot fc Falll l2000t 
IPoggianti et~ani200lt IPanuzzo et al.ll2003t) make it difli- 
cult to extrapolate the results to the galaxy as a whole 
or to other wavelengths. ICalzetti et al.l l(l99^ proposed a 
relation to link the Balmer decrement to the extinction 
in the UV and optical in starbursts, which was s ubse- 
quen t ly found not to be valid in normal galaxies ijBeli 
l2002HBuat e'ra L''2002l)^ 

V arious models (.Witt fc Gordoijr2000l: IPanuzzo et alJ 
|2003|) have shown that the FIR/UV ratio is the most reli- 
able measure of the UV attenuation, depending weakly on 
the relative geometry of the stars and the dust, the extinc- 
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tion law, or the nature of the underlying stellar popula- 
tion. This results from the fact that the stellar population 
heating the dust (emitting in the far infra-red) is the same 
po pulation responsibl e for the UV emission. 

iBuat fc Xul l)l996l) used a radiative transfer model to 
derive the extinction of 152 disk galaxies observed both 
in the UV and in the infrared, using the FIR/UV ratio. 
They found relatively moderate extinctions (0.9 and 0.2 
magnitude for early-type and late-t ype disk galaxies re- 
spectively) . It has also been used by iBoselli et alJ |2QQ^ 
in a large sample of 118 Virgo galaxies. They found an 
average UV extinction of 1.28, 0.85, 0.68 for galaxies of 
type Sa-Sbc, Sc-Scd, Sd-Im-BCD respectively. In associa- 
tion with a geometrical model, the FIR/UV ratio can be 
used to predict the dust extinct ion at all wavelengt hs, as 
was done in a simplistic way bv lBoselli et alJ l)2003|) . 

While the FIR/UV ratio has been used successfully 
in unresolved galaxies, it has never been studied in spa- 
tially resolved objects. While many profiles arc available 
for nearby galaxies (surface brightness at various wave- 
lengths, gas, metallicity), we do not know generally the 
radial distribution of extinction in normal late-type spi- 
rals. 

We propose in this paper to adopt the FIR/UV as an 
extinction diagnostic for spatially resolved nearby galax- 
ies. This is done by computing FIR and UV radial profiles, 
and deducing from them reliable extinction profiles. The 
FIR/UV ratio does not trace extinction on small scales 
(because of geometrical and transfer effects) but this prob- 
lem is avoided by considering azimuthally averaged pro- 
files and working at relatively low resolution. This radial 
variation of the dust attenuation will be compared to the 
frequently used Balmer decrement gradient. 

Our next goal will be to give to the reader an empiri- 
cal recipe for correcting for extinction UV to near-infrared 
radial profiles (at least in a statistical sense). This will be 
achieved by studying the dependence of the dust extinc- 
tion on other local properties that are likely to affect the 
amount of dust: the metallicity and the gas surface den- 
sity. 

The data used arc presented in Sect. 2: FOCA UV 
images of 6 nearby galaxies with their IRAS FIR counter- 
part. We also present the gaseous profiles and HII region 
(abundances and extinctions) data used in our investiga- 
tion. While our input data (UV and FIR fluxes, metallic- 
ity, gas surface density) tend to decrease with the radius, 
it is still unknown what is the radial variation of the dust 
extinction (as traced by the FIR/UV ratio) and of the dust 
to gas ratio. In Sect. 3, we present the extinction profiles 
obtained in the UV, and predicted at other wavelengths 
with the help of a simple model. In Sect. 4, we compare the 
attenuation in the UV with the one derived from Hydrogen 
lines in HII regions. In Sect. 5, we study the dependence of 
the extinction and the dust-to-gas ratio on the metallicity. 
We propose a simple prescription to estimate the extinc- 
tion profile in any galaxy with either a known abundance 
gradient, or a blue absolute magnitude and scale-length. 
Our most important results are summarised in Sect. 6. 



Table 1. UV characteristics of the FOCA images. 



Galaxy 


Resolution 


Exposure 




arcsec 


sec 


(1) 


(2) 


(3) 


M33 


20 


9 X 150 


M51 


12 


4 X 300 


M81 


12 


11 X 150 


MlOO 


12 


4 X 200 


MlOl 


12 


4 X 300 


M106 


20 


4 X 150 



2. Data and methodology 

The first purpose of this work is to obtain reliable dust 
extinction gradients from the FIR/UV ratio of resolved 
galaxies. This exercise is however strongly limited by the 
lack of UV images of FIR IRAS resolved galaxies, and 
by the poor FIR spatial resolution of the IRAS images. 
Indeed the sample of galaxies larg e enough to be spatially 
resolved at 100 ^im is quite small (|Ricdll99^ . 

At the present time, only six galaxies satisfy these con- 
ditions. In the UV, the images we will use were obtained 
by FO CA at 2000 A (FOCA is described in lMilhard erall 
Il99lh . Their characteristics are given in Tabled IRAS im- 
ages at 60 and 100 /im are available for 5 of them in the 
high resolution catalogue of lH,ic(j ijlOQI^ . For the last one 
(MlOO), we made an IRAS HIRES request to the IPAC 
web page^. General properties and the integrated fluxes of 
the sample galaxies are given in Table[21 These 6 galaxies 
do not constitute a complete sample in any sense, but are 
the few nearby objects for which a resolved analysis of the 
ratio FIR(IRAS)/UV(FOCA) is possible at the present 
time. 

The determination of the radial variation of the FIR 
to UV fiux ratio is done after computing the independent 
FIR (60 and 100 yum and UV radial profiles (see Sect. 
2.4). The images of each galaxy in the UV and at 100 yum 
are shown in Figs. 1 to 6 (respectively in the top-right 
and bottom-right panels). The calculation of the extinc- 
tion profiles (the UV extinction profile is shown in the 
top-left panel of Figs. 1-6) is reported in Sect. 3. 

Once an extinction gradient has been obtained from 
the FIR/UV ratio, our next goals are (a) to compare 
them with other extinction indicators (b) to see whether 
the extinction within galaxies is related to other entities 
like the gas density or the metallicity. If present, these 
relationships would be extremely useful for a more accu- 
rate computation of the effects of dust in galaxy evolution 
models in which extinction is often simply estimated from 
gas densities and abundances (e.g. Guiderdoni & Rocca- 
Volmerange, 1987). 

HII regions can provide us with additional informa- 
tion since their spectroscopic observation can be used to 
have an independent dust extinction estimate (through 
the Balmer decrement) and for the determination of the 
metallicity gradient. Data in HII regions have thus been 

^ http://irsa.ipac.caltech.edu/IRASdocs/hires_over.htnil 
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Table 2. Properties of the target galaxies. 



Galaxy 


UV 


H 


60 ^m 


100 ^m 


D 


Type 


D25 


Bt 


12 + loq(0/H) 




mag 


mag 


Jy 


Jy 


Mpc 




arcmin 


mag 


at center 


(I) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


f8) 


f9) 


do) 


M33 


6.12 


4.35 


421.6 


1102.3 


0.70 


SAcd 


70.8 


6.27 


9.08 


M51 


8.88 


5.65 


106.1 


266.8 


8.40 


SAbc 


11.2 


8.96 


9.38 


M81 


8.97 


4.09 


42.7 


161.9 


3.63 


SAab 


26.9 


7.89 


9.35 


MlOO 


10.56 


6.81 


25.9 


69.2 


17.00 


SABbc 


7.4 


10.05 


9.36 


MlOl 


7.99 


5.80 


72.1 


200.8 


7.48 


SABcd 


28.8 


8.31 


9.06 


M106 


10.02 


5.71 


26.4 


77.8 


7.98 


SABbc 


18.6 


9.10 


9.09 



Column 1: Galaxy name. Column 2: the UV magnitude from FOCA images integrated within the largest radius used for this analysis. The FOCA magnitude is 
defined as m,uY=-1.h log(i^)-21.175, where F is the flux in erg cm~ s~ A~ ). Column 3: H band magnitude IJarrett et Column 4 and 5 : Flux 

densities (IRAS) at 60 and 100 /im in Jy IHicdll9^ . Column 6 : Distance in Mpc. References: IMadore et al J 119851) ■ ^b'eldmeier et al J 11991) . Il'i:eedman et all 
rarese et al.l I200C1) . 17 Mpc is adopted for MlOO, as a member of Virgo. Column 7 : Morphological type (as found in NED^). Column 8 : Major axis 
diameter at the 25.0 mag arcsec""^ isophote in B. Column 9 : total B magnitude (columns 8 and 9 arc taken from the RC3 catalogue). Column 10: Central 
abundance (deduced from the abundance gradient extrapolated to radius — 0). 



collected for each galaxy. They are presented in the middle 
column of Figs. 1-6 and discussed in Sect.s 2.1 (extinction) 
and 2.2 (abundances). 

We finally need gas profiles. They are presented in the 
middle-left panel of figures 1-6 for each galaxies, and com- 
mented in Sect. 2.3 (the UV extinction per hydrogen atom 
is shown in the bottom- left panel of Figs 1-6). 

2.1. Extinction in till regions 

We use the A{Ha) extinctions as given by the Bahncr 
decrement in individual HII regions in each galaxy. The 
references for the data arc given in Table 13 

The Ha/H/3 intrinsic flux line ratio is relatively con- 
stant, and has a value of 2.86 in case B recombination 
llOsterbrocklll974l) . The extinction in Ha can be derived 
from the comparison of the observed ratio to this intrin- 
sic value, adopti ng a Galactic extinct ion law and a dust 
screen geometry l)Leaueux et al.lll98l|) . This observed ra- 
tio can also be affected by the underlying stellar Balmer 
absorption, which can be as strong as the emission in H/?. 
For this reason, we decided to use only the data in which 
this underlying Balmer absorption had been taken into 
acco unt, usually throug h a standard correction of 2 A 
(e.g. iMcCall et al1ll98,^ . 

Whenever possible, the extinction has been re- 
computed by comparing the observed Ha/H/? ratio to the 
theoretical one of 2.86, and removing the Galactic com- 
ponent, in an attempt to homogenize the data. The dif- 
ferences with the published values are nevertheless small 
with respect to the scatter: the largest correction is 0.5 
magnitude and the median 0.11. 

An alternative dust extinction determination in HII 
regions can by obtained by t he comparison of t he radio 
continuum to the Ha flux fe.g. lCaplan fc Deharvens-1986t 
iLeaueux et al.lll98l|) . Since such data are not available for 
all our galaxies, we are unable to make this comparison at 
this time. 

The A{Ha) extinction from the Balmer decrement for 
single HII regions is shown in the middle row-middle col- 
umn panel of Figs. ^tolHI (the number of points is smaller 



than the number of HII regions given in Table 3 because 
we kept only those with corrections for the underlying 
Balmer absorption). It shows large scatter at all radii. A 
weak gradient of decreasing extinction is nevertheless ob- 
served in all the galaxies (although hardly in M81 and 
M106), as shown by the fit (dashed line). 

For a few galaxies, independent measures of the dust 
ext inction in HII region s are available: 

IScoville et al.l l|200lh used high resolution Ha and Pa 
(Paschen-a) images of M51 to study the extinction in HII 
regions. For illustration, wc add their results to the top 
row-middle column panel of Fig. |2| Each triangle corre- 
sponds to an individual HII region, in which the extinction 
has been computed from the Pa to Ha ratio. On average, 
they find a larger extinction than the one derived from the 
Balmer decrement at th e same radius. This is partially due 
to a resolution effect, as lScoville et al.l (|200l|) estimate the 
extinction in very well defined regions in their Pa and Ha 
images. The solid line shows the extinction computed from 
the Pa to Ha ratio after the fluxes have been averaged in 
10 arcsec radial bins (only using the pixels where the sig- 
nal to noise is larger than 3 and 5 for Ha and Pa). There 
is no sign of a gradient, but the data cover a limited radial 
range (i? < 4 kpc). 

For MlOl, we present in Fig. [S] as a shaded area a 
sketch of the extinc tion variation with radius found by 
IScowen et al.l (|l992t ) in 625 HII regions identified in Ha 
and H/3 narrow-band images. The scatter in their data 
is very large (see their Fig. 3), but the trend is clearly 
detected. Inside ~ 20 kpc, they find a steeper gradient 
than the one obtained with the limited number of spec- 
troscopic studies we consider, while it is quite fiat in the 
outer galaxy (where it is however based on only a few HII 
regions). 

2.2. Abundances in HII regions 

The abundances are computed for the individual HII 
regions of each galaxy. The references for the HII re- 
gions data are given in Table El As m o st ot them 
were already collected by IZaritskv et alJ l|l994l) (and 
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Fig. 1. Profiles in M33. Top row-left column : UV attenuation radial profile obtained from the FIR to UV ratio (solid). 
The dotted curve indicates the integrated extinction within the radius R. The error bar indicates the resolution (IRAS 
resolution at 100 /xm). Middle row-left column: hydrogen density profiles (neutral: dashed, molecular: dotted, and 
totahsolid). Bottom row-left column: A{UV) to hydrogen atom column density (neutral: dashed, molecular: dotted, 
and totahsolid). Middle row-middle column : Ha extinction in HII regions (different symbols correspond to different 
data references (Table EJ, and fit (dashed line). Bottom row-middle column : oxygen abundance in HII regions and 
fit (dashed line). On the right, the UV (top) and IRAS 100 /im (bottom) surface brightness images are shown. The 
contours are separated by 1 mag arcsec"^ in both images where ii=-2.5 log(F)-|-15 {F in Jy arcsec"^) for the 100 fim 
image and /xc/y=-2.5 log(F)-21.175 {F in erg cm~^ s~^ A""'^ arcsec"^) for the 2000 A image. 



used to compute gradients), the abundances of oxy- 
gen (12-|-log(0/H) are computed from the R23 indica- 
tor ([OII]3726,3729-f [OIII]4959,5007)/H/3, using their cal- 
ibration. 

For each galaxy, these abundances exhibit a clear gra- 
dient, shown as the dashed line in the bottom row-middle 
column panel of Figs.nto|5| This fit will be used here- 
after to determine the abundance at a given radius in each 
galaxy. We implicitly make the assumption that the oxy- 
gen abundance measured in HII regions is the same as in 
all the interstellar medium at the same distance from the 
centre of the galaxy. 

In the case of MlOl, we als o show the oxygen abun- 
dance gradients deduced by iScowen et alJ l)l992() from 



their analysis of 625 HII regions (dotted lines). The dif- 
ferent curves correspond to different surface brightness 
thresholds that they apply to select the HII region s. The 
advantage of the gradients of IScowen et al.l l)l992j) is the 
large number of HII regions involved, although computed 
with limited spectral information. Contrary to other stud- 
ies, they propose a two slope gradient, but the abun- 
dances are nevertheless close to those derived from the 
spectroscopy of the limited number of HII regions of Table 

m 

With high s i gnal-t o-noise spectra of 20 HII regions, 
iKennicutt et al.l l)2003|) recently derived the electron tem- 
perature for those HII regions and determined robust 
abundances which are systematically lower (by 0.2-0.5 
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Fig. 2. Profiles in M51. Same caption as for Fig.Hlexcept: top row-middle column : Ha extinction derived from Paschen- 
a to Ha ratio, individual HII regions from lScoville et all l)200l| l ftriangles). The line is obtained by averaging the Pa 
and Ha fluxes in 10 arcsecs radial bin s, taking into account only the pixels with signal to noise larger than 5 and 3 in 
respectively the Pa and Ha images of IScoville et alJ I 20flll) . 



dex) than the abundances obtained with strong lines 
and empirical c alibra tors. Figure O shows that the 
iKennicutt et alJ l)2003(l gradient is also slightly steeper 
that from the strong-lines. 



Since we do not have other data for the rest of our 
galaxies, we use the abundance gradient as derived from 
the R23 calibration (including MlOl, for the sake of homo- 
geneity). If R23 systematically produces abundances that 
are too large, our results should be corrected accordingly, 
perhaps by reducing our abundances by a few tenths of 
dex. 



The uncertainty in the abundances is usually on the 
order of ^ 0.15 dex. This is much smaller than the differ- 
ences in the metallicity measured at various radii, so that 
the value of the gradient is relatively well defined. 



2.3. Gas profiles 

The molecular profiles are computed from the I {CO) 
profiles; references are given in Table |31 When only a 
few points are available, an exponential profile is fitted 
and used afterwards. We use a conversion factor X from 
I( CO) to H2 depen dent on the metallicity, as suggested 



i ( UUj to ti2 depen den 
bv lBoseUi etall |^02): 

logX = -1.01(12 + log{0/H)) + 29.28. 



(1) 



This correlation was found for integrated values over whole 
galaxies and is here assumed to hold for radial profiles as 
well. Note that lNakai fc Kuncl 1 199fil> found similar results 
in M51 in their analysis of the radial variation of the dust- 
to-gas ratio. The value of 12 + log{0/H) used in Eq. ^is 
the one given by the radial fit of the abundance gradient 
(see Sect. lT^ . 

References for the atomic gas profiles are given in Table 
13 Although these profiles are given at various resolutions, 
the results of our analysis should not be affected since all 
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Fig. 3. Profiles in M8i. Same caption as for Fig.^ 



Table 3. References for complementary data, Position 
Angle {PA) and inclination (incl) in degrees used for de- 
termining the profiles 



HI CO HII regions 



PA Incl 



M33 

M51 

M81 

MlOO 

MlOl 

M106 



10,11,12(21) 
10,13(10) 
14,15,16(28) 
10,17(14) 
10,18,19,20(38) 
14,21(20) 



22 
17 
-28 
153 
42 
-30 



55 
20 
59 
27 
21 
63 



References f or the HI. f ll: iDeul fc van der HulstI jl987f). [2 1: 
iRandetalJ Jl992D [3]: iRotsI lll975D. [4]: IWarmey il986D . 



[5l:IWong RlitJEloi. ffilTwevers et a^J h9H(t\. Re 



for the 
[9]: IS^ 



Corbellil i2003tl f8 |: lYoung "et alj ^ (11995 
IWong fc Blit3 ^0^^~Spectral infor- 



[5] 



motion for HII regions (for abundances and extinctions). 
The number bet ween parentheses in dicates the number of 
HII regions. f lOl: iMcCall et al.l ^ll98R^ . fi ll: iKwitter fc AUeil 
IT9S1 . f1 21 : Ivilchez et a,] J Jl 98?^) fl 31 : IPiaz et al .1 dl f 



[141 ■ [Qev fc Kennicutd (11993) [151: iGarnett fc Shieldj |1987 
_[16l:IStauffer fc Bothunl (ll984^. fl 7] : 'Shiel ds et al.l (Il99lll. [181: 
">mithl I^1975^. fl9l: iRavoigtalJ (,19821 . [20]: Ivan Zee et alJ 



liggsf) . [211: izaritskv et al.l il9" 



our profiles are smoothed to the IRAS 100 /xm resolution 
in order to compare them with the derived UV attenu- 
ation. For each galaxy, the gas profiles are shown in the 
middle row-left column panels of Figs. ^to|H| 



2.4. FIR and UV profiles 

From the IRAS and FOCA images, we compute the sur- 
face brightness profiles at 2000 A, 60 and 100 jim as fol- 
lows. After removing the stars present in the UV image 
and subtracting the sky, the resolution in all images is 
changed to that of the 100 /xm maps (100 ") by convolu- 
tion with a Gaussian filter. The pixel size of the UV image 
is changed to match IRAS (15 arcsecs per pixel). 

The fiuxcs arc azimuthally averaged using the task 
ELLIPSE within IRAF'^, keeping the centre, cUipticity 
and position angle fixed for each galaxy. We adopt the 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Fig. 4. Profiles in MfOO. Same caption as for Fig. ^ 



same position angle and inclination as used in the deter- 
mination of the HI profiles, and given in Table 13 

The UV p rofiles are co rrected for Galactic extinction 
as taken from lSchleeel et a l. (199^. 

The resulting surface brightness profiles are shown in 
Fig.CKfor the FIR, we adopt n= -2.5 log (F) + 15, F in Jy 
arcsec"^). We truncate the profiles at a radius where the 
azimuthally averaged flux surface densities become lower 
than one sigma of the sky. 

Our profile of M33 is in agreement with the one pub- 
lished by Buat et al. (1994) at the FOG A spatial resolu- 
tion, but shows less structure because our data have been 
smoothed. 

The profile of M106 shows a fiattening in the outer 
part of the UV and 60 fiux. The feature must be real 
since we are well above the sky in the UV and it is located 
in the area of a ring, clearly visible in the image. At 60 
/im, the flux at the same radius is barely more than one 
asky, and is dominated by a few large peaks, not present 
in the 100 /im image. The few outer points of this profile 
are quite uncertain. 



3. Radial extinction profiles 

The extinction profile A{UV) is determined from the 
FIR/UV ratio by combining the UV and FIR radial pro- 
files (determined as described in Sect. 2.4), using the cal- 
ibration of Buat et al. (1999), in Sect. 3.1. By combining 
these results with a simple geometrical model, we give in 
Sect. 3.2 a simple recipe to extend A(UV) to optical and 
NIR wavelengths. 

3.1. Determination of the UV extinction profile 

To compute the UV extinction A(UV), we adopt an up- 
dated version of th e fit of .Buat et al.i U 99i given m 
llglesias et aTl l)2003|) : 



A{UV) = 0.622 + 1.140 log 

\ tuv J 
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Fig. 5. Profiles in MlOl. Same caption as for Fi g. [J except: Middle row-middle column : The shaded area in- 
dicates roughly the extinction gradient derived bv IScowen et alJ lll992|) from their narro w-band imaging study. 
Bottom row-middle column : The dotted lines indicate the gradients of IScowen et af for the various values 

of the surface br i ghtne ss threshold they use with their narrow-band imaging. The long-dashed line is the gradient of 
iKennic^itt et alJ l)200,'^ measured with high signal-to-noise data (and thus avoiding the use of strong-line calibrations 
as in the other studies). 



where Fpjji (W m ^ arcsec ^) is obtained from a combi- 
nation of the 60 and 100 /im fluxes, i.e.: 

Ffir = 1.26(2.58 10^2/60 + 10^^100)10-'^ (3) 

where /go and /loo arc the IRAS surface brightnesses in 
Jy arcsec"^; and Fuv is the UV flux: 

Fuv 2000 10-^2000 [Wm-2arcsec-2] (4) 

where /2000 is the UV surface brightness (erg cm^^ s^^ 
arcsec"^). 

The IRAS FIR does not include the total dust emis- 
sion. However, the calibration of Iglesias et al. (2003), us- 
ing the ISO (Infrared Space Observatory) results of Dale 
et al. (2001), takes into account the average difference be- 
tween the FIR emission and the total dust emission in disk 
galaxies. This extrapolation is nevertheless based on data 
short-ward of 100 /xm and model SED curves and does 



not take into account the radial variation of the dust tem- 
perature. Alton et al. (1998), using ISO, found that the 
cold dust is indeed more extended than the warm dust in 
their sample of disc galaxies. We assume that the ener- 
getic balance is not much affected by this radial change, 
and implicitly use an average calibration (note that our 
60/100 ratio does not exhibit a large gradient). 

The FIR to UV ratio is believed to be a robust method 
of determining ^(C/y) fe.g. lWitt fc GordonlEoOO|) . Other 
calibrations of the determination of the UV extinction 
from this ratio have been proposed by IPanuzzo et alJ 
lj2003l) for face-on and edge-on orientations. As our galax- 
ies are quite far from being edge-on (inclinations from 20 
to 63 °), we compared our results with the A(UV) ob - 
tained for the face-on calibration of lPanuzzo et al.l l)2003|) . 
The difference in A(UV) using the two methods is less 
than 10 % except for extinctions lower than a few tens 
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Fig. 6. Profiles in M106. Same caption as for Fig. ^ 



of magnitude, as in the few outer points of MlOl (for a 
radius larger than 20 kpc, see Fig. EJ, confirming that 
the determination of A(UV) is mostly independent on the 
adopted calibration. 

An implicit assumption in using the FIR/UV ratio is 
that we take into account all the UV photons which may 
have heated the dust where it radiates the energy (after 
having eventually propagated inside the galaxy). While 
this is not the case for regions as small as HII regions, it 
is certainly true for the integrated flux of galaxies. In the 
resolved galaxies, this method can still be applied provided 
that we average the flux over large enough areas, as for the 
case of our sample where the azimuthally averaged profiles 
are determined within 100 arcsec wide annuli (between 
~ 0.5 and ~ 8 kpc in our galaxies), the IRAS 100 /im 
resolution. 

In the top-left panels of Figs.ntoEl the resulting pro- 
file of the UV extinction A{UV) is shown as a solid line 
(the dotted line indicates the extinction integrated within 
the radius i?, obtained by applying Eq.|21to the integrated 
fluxes within the ellipses of semi-major axis R). 



For all the galaxies, we observe a clear UV extinction 
gradient over a radial range several times the resolution. 
Only for M51 is a secondary peak observed, whose radius 
corresponds to its companion galaxy. Indeed the compan- 
ion suffers very strong extinction since it is not visible 
in the UV image while being prominent in FIR. The UV 
and FIR surface brightness measured in the images at the 
position of the companion give a local extinction of 3 
magnitudes. 



3.2. A simple model to derive the optical thickness 
and the extinction profile at various wavelengths 



A{UV) can be scaled to A{X) at any wavelength A once 
an extinction law (Galactic, Magellanic) and a geometry 
are assumed. To simplify the for malism and to be con- 
sistent with previous studies fe.g. iBoselli et al1l2003j) . we 
assume a Galactic extinction law and a sandwich model 
with a wavelength dependent dust-to-star scale-height ra- 
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Fig. 7. Profiles of UV and FfR surface brightness at the 
IRAS 100 iiva resolution (~ 100 arcsec). For the infrared, 
^ is defined as -2.5 log {F) + 15, F in Jy arcsec"^. In 
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tio ^. The extinction A (in magnitude) then depends on 
the optical thickness (t) via Eq. El 



-2.blog 



1-e 



(1 



(l-e- 



.(5) 



Ti, Ai, and ^ are wavelength-dependent, and Ai corre- 
spond to a given inclination, i.e. = rosec(i). The relation 
between Ai and is shown in Fig. |S1 for several values of 

We numerically invert Eq. |S1 to derive the profile 
of Tuv from the one of A{UV). Then, we compute 
the optical thickness at a given wavelength as : ta = 
Tuv k{UV), where fc is a typical Galactic extinc- 

tion curve (Rv = Ay I E{B - 1/)=3.1). The choice of this 
extinction curve is justified by the fact that the metallic- 
ity in our galaxies (see Table |21 and Fig. ^1 for instance) 
is larger than in the Magellanic Clouds. We are implic- 
itly assuming that the albedo does not depend strongly 
on the wavelength and that the scattering is isotropic (we 
also tried to i nclude the effects of the albedo and phase 
function as in ICalzetti et al.l l)l994ll with very similar re- 
sults). 
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Fig. 8. Relation between the absorption A,; and the op- 
tical thickness for sandwich models of various dust to 
star scale- height ratios ^. 

We can then compute the extinction at any wavelength 
with the help of Eq. |S1 We show the result of this oper- 
ation for the V and H band wavelengths in Fig. |^ We 
adopt a dependence of the dus t-to-star scale-heigh t ratio 
on wavelength as suggested bv lBoselli et all ()2003|) : 



= 1.0867 - 5.501 10"^A(1), 



(6) 



This assumption is used to mimic the fact that the young 
stars (predominantly emitting at short wavelengths) lie in 
a thinner disk (similar to the dust, ^=1 at 2000 A) than 
older stars (emitting at long wavelengths) which migrate 
to larger heig hts with age. This is in ag reement with some 
observations l)Boselli fc Gavazzil Il994|) . b ut not always 
confirmed for edge-on nearby galaxies (e.g. IXilouris et alJ 
EM). 

Note that a sandwich model with a constant dust-to- 
star scale-height ratio ^(A)=0.5 is not consistent with the 
observations. It can be seen in Fig. |S1 that for this value of 
^, the absorption in magnitudes can reach only moderate 
values. Even for infinite optical thickness, we would ob- 
tain only A{UV)=1.5 which is lower than the absorption 
measured in the centres of M51, M81 and Ml 00. 

The results of this model applied to our galaxies are 
shown in Fig.|51for the wavelengths of the V and H bands. 

The variation of extinction with radius that we can 
detect is between 0.5 and ~ 1 mag over the whole radial 
range in the UV. The results of the model indicate that 
the typical V band extinction gradient is between 0.2 and 
0.5 mag, and it is very small in the H band. Of course this 
is model-dependent and could easily vary by a few tenths 
of magnitude in the V band if we change ^. 

Fig. 1^1 also shows the extinction profiles in V and H 
when a simple screen of Milky Way- l ike dust i s ado pted 
(marked "S" in the figure) or with the lCalzettil l)l999() law 
( "C" ) . The differences between these models are negligible 
at long wavelengths but can reach 0.3 magnitude in V (and 
eventually change with the radius). 
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This shows that the geometry (and the choice of one 
geometrical and dust model rather than another) can have 
an important impact on the determination of the colour 
gradient of the underlying stellar population, even if the 
observed colour gradient is unlikely to be mainly due to 
the dust ijde Jong|ll996|) . The extinction dependence on 
the geometry is a well known effect, studied in detail in 
Disney et al. (1989). 

3.3. Comparison with other studies of M51 

iHill et all l(l997f) show a UV colour profile of M51 (//151- 
Hu) varying by 1 mag or so between ~ 2 and 12 kpc (the 
gradient is even stronger if smaller radii are included, but 
we cannot probe such small scales with our resolution). 
Hill et al. propose that the colour gradient is due to an 
extinction gradient. We indeed find (in a more direct way) 
an extinction gradient in the UV in M51 (Fig.[51). 

IVanseviciuj ()200l|) estimates the extinction at 11 kpc 
from the centre of M51 using its companion as a tracer. He 
finds Ay = 0.7 ± 0.3 mag. The extinction we derive from 
our model at this radius is a bit lower (~ 0.5 mag) but 
within the error bars (see Fig. |5J). As another illustration 
of the dependence of our results on the adopted geome- 
try, changing ^ to a constant value of 0.7 (independent of 
the wavelength) would increase the value of Ay by only 
0.1 magnitude (larger differences would occur where the 
optical thickness is larger, i.e. in the inner few kpc of the 
galaxy). 

4. Extinction in the UV continuum and in HI! 
lines: A(Ha) vs A(UV) 

To check the consistency of our extinction determination 
with average values available in the literature, we estimate 
the total integrated extinction for both the UV and Ha 
and compare them to the results of lBuat et al.l 1 20021) for 
late- type normal spirals. 

The integrated UV extinction within the radius where 
it is computed can be seen as the dotted curve in the top- 
left panels of Figs. to El It usually presents a plateau 
or a small decrease at large radii; we then consider as 
"integrated" extinction over the whole galaxy its value at 
the largest radius. We choose as a typical Ha extinction 
the average of A(Ha) found in the HII regions considered 
in this study. 

The average Ha extinction (and its standard devia- 
tion) and the UV integrated extinction are compared in 
Fig. ^1 This figure also shows the local values (at differ- 
ent radii) within each galaxy connected by a curve. The 
integrated value is representative of an "average" of the 
local values obtained at differe nt radii. Th e filled circles 
are the star-forming galaxies of iBuat et alJ ll2002l) . using 
the updated Ha extinction of lGavazzietaLl <)2004|) (which 
introduces small differences with respect to Buat et al.). 

Our "average" values (squares) and most of the profile s 
(curves) largely overlap the points of iBuat et al.l ((20021) . 
however, we span a limited range of extinction than they 



do since we are limited to absorptions less than 2 mag, 
while Buat et al. include objects with up to ~ 4 magni- 
tudes. Actually, this is the case of MlOO for which Buat et 
al. find A{Ha) > 3.7 mag. However, they notice that this 
high value may be due to the very high surface brightness 
nucleus while the disk might be less extincted. Our Ha 
extinction derived from HII regions in the disk is indeed 
lower. 

An important d ifferen ce between our approach and 
that of iBuat et al.l l)2002|) is that we defined an average 
A{Ha) in HII regions, while their integrated value is dom- 
inated by high surface brightness regions and the nucleus. 
This might explain the differences between the two studies 
and especially why o ur "local" study do es not reach the 
high values found by iBuat et'all l|2002l) in "integrated" 
galaxies. 

We can also see in Fig. ^] that lower values of the 
integrated UV extinction are found for M33 and MlOl 
(which are of type Scd) than for the other galaxies (type 
Sbc and Sab). The average values for these two groups are 
in agreement with (albeit slightly l ower than) the average 
values given bv lBoselh et"an j2003(l for Sc-Scd (0.85 mag), 
and for Sa-Sbc galaxies (1.28 mag), obtained for a much 
larger number of integrated galaxies. 

The profile of M81 in A{Ha) is rather flat, as seen 
in Fig. 1101 It can be seen in Fig. 13 {middle row-middle 
column panel) that this is due to a large scatter of the 
extinction in a limited number of HII regions. In general, 
the average values of A{Ha) (squares in Fig. llO|l arc more 
robust than the local values as the latter are more sensitive 
to the small number of HII regions. 

Finally, we also show in Fig. Uni the ICalzettil lll997^ 
att enuation law for starbursts {A{UV) = 1.6A{Ha), 
see iBuat et all I2OO2L for the derivation) The relation is 
st eeper for the integ rated star-forming galaxies, as shown 
in lBuat et al.l l|2002|) . although consistent, within the large 
uncertainties, with local values. 

5. Extinction, metals, and gas 

5.1. Dependence of the extinction on the metallicity 

iQiiiUen fc Yukital 1I2OOJ) . using extinction derived from 
the Paschen-a (Pa) to Ha ratio in the centres of a few 
galaxies and along the profiles of M51 a nd MlOl (with 
original data from the thesis of IScowent .199211 , found 
that the extinction correlates with the metallicity and 
with the gas surface density (their Fig. 4). lHeckman et all 
(il998) have shown a dependence of the UV spectral slope 
P on 0/H in their starbursts galaxies (0 is defined as 
F(X) (X for 1250 < A < 2600 A by iCaJzetti et a,l1 
(|l994() and is equal to ^ -2.1 for a dust-free starburst). 
Several studies, theoretical and observational, have tried 
to relate the amount o f extinction to m etallicity and gas 



Issa et al. 1990- .Liscnfcld fc Ferraralll998l: 



amount (e ^ „ .. 

lDweklll99Sl:lHira,shita et a,1.ll200l|> 



From eq. El (adopting the varying dust-to-star scale- 
height ^, as given in Eq. EJ we derive for each profile of 
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Fig. 9. Extinction profiles in the UV, V and H. A{UV) is obtained from the FIR/UV ratio. The black Une shows the 
extinction profiles in V and H as derived for a sandwich model with dust to star scale-height ratio ^ depending on 
the wavelength. For each of thes e two bands, th e others curve show the predicted extinction for a dust screen with 
Milky- Way- type dust (S) and the ICalzettJ (|l99i3l) law (C). 



A{UV) the corresponding profile of Tuv.i, which is cor- 
rected for inclination, to obtain rj/y^o- 

In Fig. (bottom), we superpose the relation be- 
tween the optical thickness tuyq and the abundance for 
all the galaxies. In the top panel, we directly plot the 
extinction A{UV) in magnitudes (uncorrected for inclina- 
tion) as a function of the oxygen abundance. The advan- 
tage of using A{UV) as derived from the FIR to UV flux 
ratio with respect to other dust extinction indicators is 
that it does not depend on the extinction model. A quite 
good correlation between the extinction and the metal- 
licity is observed, suggesting that the extinction does de- 
pend on the metallicity. This dependence is confirmed by 
the U V extinctions of th e integrated star forming galax- 
ies of iBuat et al shown as circles, adopting the 
metallicities of Icavazzi et al.l (|2004i and in preparation), 
although with a larger scatter. This correlation between 
local extinctions and metallicities amongst our different 
galaxies is a striking result as the oxygen abundance and 
A{UV) are obtained independently of each other. 



Actually, a dependence of the UV extinction on 
themetallicity was already suggested bv iHeckman et all 
()l998t) for their starburst galaxies. We show their data a s 
triangles in the same figure. Following iBuat et al.l 1 20021) . 
the A{UV) extinction at 2000 A in starbursts is computed 



A{UV) = 0.9 X 2.5 log I 



0.9 F , 



(7) 



UV,1600A 



The extinction in starbursts is much larger for the same 
metallicities and shows a larger scatter than what we find 
within the disks of our spirals, as already remarked by 
iRuat et a1.l|200i. 

The black lines are the median of the fit performed for 
each of our galaxies. The numbers given in parentheses in 
eq. 8 indicate the standard deviation in the parameters 
found for each galaxy. These numbers reflect the large 
differences between galaxies, nevertheless the global trend 
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1 2 3 

A(UV) (mag) 

Fig. 10. UV vs Ha extinction. The squares are the av- 
erage values for each galaxy. The curves are local values 
at different radii within each galaxy . Filled circles are the 
star-formii ig galaxies of iBuat et alJ l(2002l ) (adopting the 
A{Ha) of iGavazzi et al.ll2004 - computing A{UV) from 
the UV/FIR ratio, like for our galaxies, but using inte- 
grated fluxes. The long dashed line indicates the relation 
A{UV) = 1.6A{Ha), expected in the case of the attenua- 
tion law proposed bv .Calzettii ((l997l) for starbursts. 



is well established. 

logiTuv.o) -15.25(±6.27)-H.78(±0.70) x Z 

or : (8) 
A{UV) = -21.80(±11.70)4-2.54(±1.28) X Z 

where Z = 12 + log{0/H). The standard deviations of r 
and A{UV) around the median-fit are cr=0.37 and 0.43 
respectively. 

The large scatter produces an important uncertainty 
in the zero-point of this relation, even though the depen- 
dence on the metallicity is well defined. These relations 
could be used to predict r (or A{UV)) as a function of the 
radius for each galaxy with an available abundance gradi- 
ent. In the case where no abundance gradient is available, 
we can use the fact that the abund ance gradient is con- 
stant in units of the disc sca le-length ijPrantzos fc Boissieil 
l20nnHHenrv fc Worthe\ll99ft. and references within). The 
oxygen gradient is typically -0.2 dex / Rb {Rb is the disk 
scale-length in the blue band). The zero point of the re- 
lation can b e determined fro m the magnitude-metallicity 
relationship. iGarnetd l)2002|) has shown that the oxygen 
abundance at the effective radius in spiral galaxies satis- 
fies: 




8.5 9 
12 + log(0/H) 



9.5 



Fig. 11. UV Extinction (top) and optical thickness (bot- 
tom) as a function of 12+log(0/H). Each galaxy is rep- 
resented by a curve. The symbols correspond to radii 
separated from each other in units of the spatial reso- 
lution. The median of a fit performed on each galaxy 
is shown as black line. In the top panel, we also show 
the extinction derive d for the integrated starbursts of 
iHeckman et al.l (Il998l) as tri angles and the star forming 
galaxies of lBuat eT^ l)2002() as circles. 



where AIb is the blue absolute magnitude. For an expo- 
nential disk, Reff = 1.865i?B- Combining this result with 
Eq. |S1 we obtain a second expression of the optical thick- 
ness with no explicit dependence on the abundance gradi- 
ent: 



log{Tuv,o) 

or : 

A{UV) 



= -0.28(±0.11)Afi3-0.36(±0.14) 
- 4.62(±10.45) 



R 

Rb 



(10) 



0.40(±0.20)Mb - 0.51(±0.25) 
6.60(±19.38) 



R 
R^ 



[log{0/H)]Reff = -OAGMb - 6.4 



(9) 



These relations can be used for any galaxy with an 
absolute magnitude Mb and scale- length Rb- The uncer- 
tainties given in Eq. 1101 are obtained only by propagating 
the uncertainties of Eq. 
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Fig. 12. Top: Radial profiles of the difference between 
A{UV) deduced from the FIR/UV ratio and A(UV) 
obtained by applying Eq. 8 to the oxygen abundance. 
Bottom: radial profile of the ratio of the optical thickness 
deduced from the FIR/UV ratio to the one obtained by 
applying Eq. 8 to the oxygen abundance. The horizontal 
lines show the average value and ± one tr. 



To check the validity of Eq. 8 for computing an extinc- 
tion profile from the abundance gradient, in the top-panel 
of figure 12, we show as a function of radius the difference 
between the extinction deduced from the FIR/UV ratio 
(^(C/F)i?/ij/(7y) and the one deduced from the oxygen 
abundance at this radius (A([/U)o///). The bottom panel 
shows the ratio of the two optical thicknesses. 

The optical thickness "predicted" by applying Eq. 8 
to the oxygen abundance, tq o/h^ is within a factor two 
of the "observed" one {t^ pij^^uy) for most of the points. 
Large departures are found only for the case of M33, which 
is the galaxy with the smallest metallicity in our sample. 

Similarly, the A{UV)o/h deduced from the metallicity 
is frequently found at ~ 0.3 mag from A{UV)fib,/uv ^-nd 
the largest deviation is for MSB (up to more than one 
magnitude in its outer parts). This may indicate that our 
relations are valid only in the high-metallicity domain. 

The points situated at radii larger than ~ 20 kpc also 
show large differences between the predictions and the ob- 
servations. However, these data only concern the outer 



part of one galaxy (MlOl) where A{UV) is very small 
anyway (and where the calibration of our method is not 
as secure: see Sect. I3.l(l . In any case, deviations from the 
average relation in the external parts are expected since 
the uncertainties on the profiles (especially in the FIR) 
are larger. 

5.2. Dust-to-Gas ratio 

The extinction-to-gas ratios [bottom-left panel of Figs. ^ 
to|ni) ^'''6 computed from the gas profiles [middle row-left 
column), after smoothing the resolution to the IRAS 100 
//m maps. This ratio shows various values and trends for 
each galaxy. 

Based on extincti o n me asured in HII regions at var- 
ious radii, llssa et all (|l990() suggested that the dust-to- 
gas ratio depends on the metallicity. This is quite natural 
since the dust is formed from metals and the mass of met- 
als is equal to the product of the metallicity and the gas 
mass. This trend is naturally obtained in the framework 
of models computing consistently the evolution of metals 
and dust, despite the l arge uncert ainties in the yields of 
bo th (e.g.lT^ouell2003trnwd3ll99St). 

iGuiderdoni fc Rocca-Volmerangj l)l987(l used the the 
solar neighborhood and the Magellanic Clouds to propose 
the relation tx/Nh oc [Z/ZqY with s - 1.35 (1.6) for 
wavelengths shorter (longer) than 2000 A, [N}j represents 
the total hydrogen column density). 

As can be seen in Fig. ^| (top), we obtain a similar 
trend, although with a large dispersion, especially for the 
mctal-rich part. A formal fit (short dashed line in Fig. 
I13|l to our data gives a slope of 0.88. This is lower than 
the va lues proposed bv iGiiiderdoni fc R,occa-Volmerang3 
l|l987t) . which are illustrated by the two long-dashed lines 
in Fig. ^1 Note however that we apply different methods 
(FIR/UV ratio, metallicity dependent CO to H2 conver- 
sion factor...) and probe different ranges of metallicities 
as our spirals are more metal-rich than the Magellanic 
Clouds. 

Following iHirashita et al ] (l200j) . we can convert 
our optical-thickness-to-hydrogen ratio into a dust-to-gas 
mass ratio. This is done by using the profiles of Ay (Fig. 
El computed with the model of Sect. 13.2(1 and assuming 
that Md/Mq (the dust-to-gas mass rati o) is equal to 6 
xK)-3 when Nh/E[B - V)^5.9 lO^^ jHirashita et alJ 
l200l[ and references within) . The results are shown in the 
middle panel of Fig. El In this figure, we comp are our re- 
sults w ith the dust-to-gas mass ratio found bv llssa et alJ 
(119901) for 6 galaxies. They were deduced from the Ay/Nn 
of HII regions (or galaxy counts and interstellar absorp- 
tion in the Milky Way), estimated at a r adius of 0.7 Rdv 
[Rdv is from Ide Vaiicouleu rs et^l 1 97fil) . 

Three of the galaxies of llssa et al.l (|l990|) are in com- 
mon with our study (M33,M51,M101), presenting some 
differences with our results that must come from the fact 
that we derive Ay from the FIR/UV ratio and use the 
gas profile, while they use measurements done in HII re- 
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Fig. 13. Top: Optical thickness to hydrogen column 
density ratio. The short-dashed line is a fit to our 
data of slope 0.88. The long-dashed lines show the 
metallicity dependence of the dust to gas ratio of 
iGuiderdoni & R.occa-Volmerangel l(l987f) . middle: dust to 
gas mass ratio (computed from the Ay/Nn ratio as ex- 
plained in Sect. 15. 2|) as a function of metallic i ty. Th e cir- 
cles mark the integrated galaxies o fllssa et al.l (llQQOl). Th e 
solid curve is the best-fit model of lHirashita et al.l l(200ll) . 
bottom: dust to gas mass ratio (computed from the 100 /im 
fl ux density). The cr osses indicate the integrat ed galaxies 
of lTuffs et aLl l|2002l) and lContursi et al.N200l[l. The solid 
curve is the best-fit model of lHirashitgretaLT i 200l|) . 



gions. In MlOl, our metallicity is higher, but our gradient 
is based on several, more recent studies. 

A more direct way to estimate the dust-to-gas mass ra- 
tio is to deduce the dust mass from the FIR flux. We com- 
pute the dust mass following iDevereux fc You nS ill 9901) 
with the same numerical coefficient as in iBosefliet'ldl 
l|2nn2l) : 



Md = 1.27Fiooi?'(ea;p(144/rz5) - 1)Mq, 



(11) 



The mass is very sensitive to the dust temperature 
which is still poorly constrained by the 100 jum IRAS 
data. In the bottom-panel of Fig. we show as crosses 
the dust-to-gas mass ratio obtained from the ISO measure- 
ments of integrated flux densities of galaxies bv pTuffs et all 



(l200jl^ andlContursi et ap l)200l|) , where the metallicity is 
from lGavazzi et alJ l)2004() . We ado pt the average tem per- 
ature (Tr,=18.9 K) of the sample o flTuffs et aLl (|2002tl (for 
which the dust temperatures are given in Popescu et al., 
2002) to compute the dust-to-gas mass ratio of our galax- 
ies at different radii. This temperature is consistent with 
the 200 ^m observations of Alton et al. (1998). The result 
is shown in the same figure as a function of the metallicity. 
The shaded area indicates a dispersion of ± 3.7 K. 

The dust-to-gas ratio is lower than the one deduced 
from the Av/Nh and shows a clearer trend with metal- 
licity, still with a large dispersion. 

If a dust temperature gradient is present (with higher 
To in the center of galaxies than in their outskirts) , as ex- 
pected from the observed gradient in metallicities and star 
formation rate (and suggested by the ISO observations of 
e.g. Alton et al., 1998), it would flatten the relation be- 
tween the dust-to-gas mass ratio and the metallicity, mak- 
ing it more similar to the one derived from the Ay/Nn 
ratio (see bottom panel). 

The best-fit model of iHirashita et al is indi- 

cated in the figure by the solid curve. It reproduces the 
global trend of increasing dust-to-ga s ratio from dw arfs 
(|Lisenfeld fc Ferraralll998|i to spirals lllssa et al lll990l). 

W hen a dust-to-gas ratio is determined as in llssa et ali 
lll99(t . the trend for our galaxies is weaker than in this 
model, and a great deal of diversity is observed amongst 
the galaxies. The values obtained from the 100 /xm surface 
brightness are in slightly better agreement with the model. 
The uncertainties due to the temperature (shaded area) 
and its gradient (see discussion above) are however large. 

6. Conclusion 

We have combined 2000 A UV images obtained with 
FOCA, and FIR IRAS images at 60 and 100 ^m to com- 
pute the FIR and UV profiles of six nearby late-type galax- 
ies. We used the FIR/UV ratio to trace the radial variation 
of the UV extinction in each galaxy. We detect a mono- 
tonic gradient of decreasing extinction with the radius in 
all of them, except in M51 where the companion produces 
a second peak in the averaged profile. 

These extinction profiles were compared to the ex- 
tinctions derived from hydrogen lines (mainly the Balmer 
decrement) in HII regions and we studied the relation 
between the extinction, the gas surface density, and the 
metallicity. 

The most significant result of this analysis is a clear 
correlation between the UV extinction (in magnitudes or 
in optical thickness) and the metallicity deduced indepen- 
dently from the FIR/UV profile and the abundance gradi- 
ent respectively. This correlati on is also found in the inte- 
grated star forming galaxies of iBuat et alJ ll2()( ^2| ) fa sim- 
ilar r elationship was found in starbursts (|Heckman et alJ 
Il998|) . but with larger scatter and higher extinction for a 
given metallicity). 

A fit to our azimuthally averaged data provides a sim- 
ple relationship between the extinction in the UV and the 
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metallicity in galaxies for which an abundance gradient 
is available (Eq. ISJ. Coupling this result with the mass- 
mctallicity relationship, wc derived an empirical formula 
linking the extinction profile to the blue absolute magni- 
tude and disc scale-length of galaxies (Eq. ^J). Once the 
UV extinction profile is determined by one of the previous 
methods, the extinction profile at any wavelength can eas- 
ily be derived through a simple model as the one described 
in Scct.O 

Our results were obtained with a very limited sample 
of galaxies. In the future, we will however be able to extend 
this work with the new images of GALEX in the UV and 
SIRTF in the far-infrared for a larger number of galaxies. 
GALEX data are being used to study the extinction radial 
profile in individual spirals: in M83 (Boissier et al., 2004) 
and in a comparative study with ISO FIR maps for MlOl 
(Popescu et al., 2004). 
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